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Abstract. Wilson's Numerical Renormalization Group (NRG) is so far the only 
nonperturbative technique that can reliably access low-energy properties of quan- 
tum impurity systems. We present a recent extension of the method, the DM-NRG, 
which yields highly accurate results for dynamical quantities at arbitrary frequen- 
cies and temperatures. As an application, we determine the spectrum of a quantum 
dot in an external magnetic field. Furthermore, we discuss magnetic impurities 
with orbital degeneracy, which have been inferred in recent experiments on quan- 
tum dots in an Aharonov-Bohm geometry. It is demonstrated that for spinless 
electrons, interference between neighbouring levels sets the low-energy scale of the 
system. Switching on an external field leads to a remarkable crossover into a regime 
dominated by orbital Kondo screening. We predict that the broadening-induced 
level splitting should be clearly visible in measurements of the optical absorption 
power. A more general model including the electron spin is studied within an ex- 
tended two-band NRG procedure. We observe competition between interference 
and Kondo screening, similar to the situation in two-impurity models (RKKY). 



1 Introduction 

Quantum impurity models and their low-temperature properties are of cen- 
tral importance in condensed matter physics. They show characteristic many- 
body effects like the screening of a local moment by conduction electrons (the 
Kondo effect) which was first observed in measurements on dilute magnetic 
impurities in metals (see ]l|). More recently, artificial nanostructures (quan- 
tum dots (^] or surface atoms probed by STM ||Q) with tunable param- 
eters provided new representations of the Anderson or Kondo model fi@|. 
In theory, a very fruitful line of research was opened by the development of 
dynamical mean-field theory (DMFT) |?J where correlated lattice systems 
are mapped onto effective impurity models which are then accessible in a 
controlled way || . 

In this article, we focus on semiconductor nanostructures, where - at the 
moment - experiments with the highest level of control can be performed. 
Electronic transport through ultra-small quantum dots, where the charging 
energy is the largest energy scale, has been studied extensively over the last 
few years Due to the quantization of charge the transport is dominated 
by Coulomb blockade. More recently, experiments revealed that the Kondo 
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effect leads to an enhancement of the conductance - the zero bias anomaly — 
in the Coulomb blockade regime dJl(| as predicted some time ago ADJ. 

Theoretical modelling of these systems is usually based on the Kondo 
or Anderson Hamiltonians |^],[| describing a localized spin (orbital) which 
is coupled to one or several conduction electron reservoirs. In the regime of 
interest, this coupling is the usually the smallest energy scale. It was realized 
very early |6| that a treatment of these models based on perturbation theory 
fails due to logarithmic divergences below a characteristic temperature scale, 
the Kondo temperature Tk- 

Solution of the Kondo problem at T <C Tk thus required a non-perturbative 
technique, which was provided by Wilson's pathbreaking Numerical Renor- 
malization Group (NRG) |p"2f . It proved to be very successful in clarifying 
the low-energy properties of various impurity problems |^3| , pT|Jl^ , [l^ , |l7t , and 
it will be the method of choice for analyzing more complex quantum dot 
systems. 

In this article we give a short introduction to the NRG technique, includ- 
ing in particular a recently developed density-matrix formalism (the DM— 
NRG |Q) necessary for a reliable calculation of dynamical quantities. Using 
the new algorithm, impurity spectra are calculated for a quantum dot in 
an external magnetic field. Finally, we discuss the subtle interplay between 
interference and interaction which arises in a quantum dot with orbital de- 
generacy. 



2 Generalized Numerical Renormalization Group 

In the following, we consider the Anderson Hamiltonian |^j 

#and = £k ^ C k<x+E Vk if* C kcr + h - C ) + U n /T n fl+ £ f n f~ h Sf(l) 

kcr kcr 

where the hybridization T(u)) = 27r^ fc |Vfc| 2 5(lu — tk) between the / impurity 
and the reservoir electrons Cko- is balanced by a local Coulomb repulsion 
U which suppresses double occupancy of the impurity. In addition, a local 
magnetic field h is coupled to the impurity spin Sj. Units are chosen as 
h = ks = /is = g = 1 and the half bandwidth is given by D = 1. 

The key idea introduced by Wilson is the logarithmic discretization of the 
conduction band shown in Fig. |J where each energy scale is represented by 
a single fermionic degree of freedom. After performing a Lanczos transforma- 
tion (for details see p^ , ^9| ) the conduction band can be written as a linear 
chain shown in Fig. 




d Ld {n+1)a + h.c.) (2) 



with hopping coefficients decaying exponentially as e n ~ A n l 2 . In this rep- 
resentation, the impurity is only coupled to the maximally localized reservoir 
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state do = (1/yN) ^2 k c^. The transformed model - while still a nontriv- 
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Fig. 1. Logarithmic discretization of the conduction band. 
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Fig. 2. Linear chain representation of the Anderson impurity Hamiltonian. 



ial many-body problem - can now be solved by iterative diagonalization, 
keeping in each step only the lowest, most relevant levels. This procedure 
resembles the one employed in calculating atomic spectra and is illustrated 
in Fig. |J Additional symmetries like the conservation of the total charge and 
components of the total spin can be invoked in order to simplify the remain- 
ing matrix algebra. In the single-band case considered here, matrix size is 
not a problem, while for calculations with two (different) reservoirs the use 
of symmetries can be vital in order to render the problem manageable. 
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Fig. 3. Iterative diagonalization of the impurity problem. 



In the iterative diagonalization scheme, the number of iterations corre- 
sponds to the temperature one is interested in according to Tn — c J 4 _ ( Ar_1 )/ 2 ) 
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where c is a constant of order one. For calculating static thermodynamic ex- 
pectation values, all necessary information is thus obtained because only ex- 
citations on the scale TV are relevant. As an example, consider the impurity 
magnetization 

< Sf > T = Z- 1 e ~ PE " < m\Sf\m > N (3) 

m 

where the \n > are the many-particle eigenstates of H and Z is the partition 
function. Due to the Boltzmann factor, higher excitations - already lost in 
iteration N - can be safely neglected at this point. 

The situation changes completely when we consider a dynamical quantity 
like the spin-resolved spectral density 

^ e -0E m , -0E„ 

A <7 (u J )=J2\<m\fl\n>\ 2 S(L0-E m +E n ) ± . (4) 

nrn 

Obviously, spectral information at frequencies u> ^> T/v requires matrix el- 
ements between low-lying states and excitations which in iteration N are 
not available anymore (they have already been lost by truncation). In or- 
der to deal with this situation, the following two-stage procedure has to be 
employed: 

(1) NRG iterations are performed down to the temperature Tjv of interest, 
in particular we choose TV -C Tk to calculate ground-state properties. In 
each iteration step, we keep the information on the transformation between 
one set of eigenstates and the next, i.e. we save the corresponding unitary 
matrix. After obtaining the relevant excitations at temperature TV one can 
define the density matrix 

p= Z _1 ^e _ - B ™/ T " I to > N < to I (5) 

m 

which completely describes the physical state of the system. The equilibrium 
Green's function can be written as 

G T (t)=^(t)Tr(p{/ T (t),/)(0)}) (6) 

(2) Now we repeat the iterative diagonalization for the same parameters. 
Each iteration step N' yields the single-particle excitations (and matrix el- 
ements of p) relevant at a frequency oj ~ Tjy/. But instead of using (Q), 
we now employ (0) and evaluate the spectral function with respect to the 
correct reduced density matrix |2(| : As depicted in Fig. |J, the complete chain 
is split into a smaller cluster of length N' and an environment containing the 
remaining degrees of freedom. In the product basis of these two subsystems, 
the full density matrix has the form 

P= Pminumznzlmi >env |^l > sys < n 2 | < m 2 | (7) 

mim 2 
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< P > 

Fig. 4. Reduced density matrix obtained by tracing out "environment" degrees of 
freedom of the chain. 

which is in general not diagonal anymore. Performing a partial trace on the 
environment then yields the density submatrix 

^red = p" d „ 2 |m > sys < n 2 | , p^ d „ 2 = Pmnumn, (8) 

nin2 m 

This projection is easily done using the previously stored unitary transforma- 
tion matrices. Note that p red - defined only on the shorter chain - contains 
all the relevant information about the quantum mechanical state of the full 
system. 

The single-particle spectrum calculated in this way is shown in Fig. ^ 
With increasing magnetic field, the Kondo resonance is suppressed and even- 
tually merges with the lower atomic level. Regarding the total density of 
states (DOS) A(ui) = Ao-(u>), the Kondo peak is split by the field and 
the DOS at the Fermi level strongly reduced. This effect has been observed 
directly in measurements of the differential conductance through a quantum 
dot§. 
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Fig. 5. Spin-dependent (left) and total (right) impurity spectral density at zero 
temperature for r — 0.02, U = 5F, and e/ = — 2.5F. The Kondo temperature is 
T K = 6.8 x 10 -4 . 
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3 Interference and Interaction in Multi-Level Dots 

After the discussion of the spin-degenerate Anderson impurity model ([j]) in 
the last section, we now consider the effect of orbital degeneracy. We will first 
study a dot consisting of two levels without spin or, equivalently, two dots in 
an Aharonov-Bohm (AB) geometry with one level per dot in the presence of 
an interdot Coulomb repulsion U. Such a system is of fundamental interest 
since the two possible paths through the dot (via level 1 or 2) can interfere 
with each other. The interference can be controlled by an AB flux and has 
attracted much interest due to the possibility of realizing AB interferometers 
[ 22 or using the coherent properties in connection with quantum computing 
[23 . Furthermore, there is enhanced experimental interest to study quantum 
dots in the strong tunneling regime where the level broadening is of the order 
of the level spacing. In this case, transport is inevitably controlled by multi- 
level physics. 




Fig. 6. Two quantum dots in an Aharonov-Bohm geometry. 



Let us first discuss the case of spinless electrons, assuming e.g. a large 
Zeeman splitting p4j . The Hamiltonian is then written as 

H = J2e kr al r a kr + ^(^ f 4r c J + h - c + ^Z £ J c ] c i + Un l n 2 (9) 

kr rkj j 

where j = 1, 2 labels the two levels and the dot is connected to two reservoirs 
r = L, R via tunnel barriers. Note that the index labelling the dots is not 
present in the reservoirs - this model contains no conserved quantum number 
corresponding to spin, unlike previous studies p5| , p6|j2^ , p8| . The tunnel ma- 
trix elements are assumed to be real except for an AB-phase, i.e. we attach a 
phase factor e 1 ^ to . The level broadening is defined by IJ — 2n\Vj\ 2 po, 
where po is the density of states in the leads, which we assume to be constant 
in the energy range of interest. The total broadening of each level is therefore 
given by r = Ff + rf. 

Since both levels overlap with the reservoirs, they have an effective overlap 
matrix element A, which induces a level splitting Se = \J 5e 2 + \ A\ 2 where 
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Se = €2 — ei denotes the bare level spacing. Within second-order perturbation 
theory it is established that A vanishes for the noninteracting dot, while 

in the case of strong on-site repulsion U ^> \e\, r. For positive level energies, 
the tunnel splitting can be observed directly as a shift of the upper level 
position (see Fig. 0) in the total spectral density 

A(u) = -- y ImG«(w+). (11) 

For low lying levels e < —T and a large Coulomb repulsion, the dot is singly 
occupied at low temperatures. In this case, the effective level splitting shows 
up in a many-body resonance ( "shoulder" ) in the spectral density at a pos- 
itive frequency u> ~ Si, see Fig. This new many-body energy scale can 
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Fig. 7. Left: tunnel splitting visible in the total spectral function for parameters 
U = 16-T, ei = <2> = T = and different positions of the second levels £2- Right: 
Absorption power (broken lines) vs. single-particle spectral density (full lines') for 
U = 50r, ei = -10r and * = T = 0. In the inset, U = 10r and ei = -3.5r Q. 



be seen most clearly in a microwave absorption experiment, where tran- 
sitions between the two dot levels are induced due to the dipole operator 
O = c i c 2 + c 2 c i- The corresponding spectral density 

p aha {u) = Z-'Y,\ <n\0\m> \ 2 S (lo + E n - E m ) (e"^» - e^-) (12) 
win 

displays a well-pronounced resonance at the frequency <5e (see Fig. ^). 

For (/> = 7r and J 1 ^ = J^, the tunnel splitting is zero, and the system is 
shown to be equivalent to an Anderson model with Zeeman splitting Se. This 
can be seen most easily by introducing the new levels 

(l/v/2)(ci±C2) (13) 
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which are then coupled to the right and left reservoirs, respectively. In this 
way, an orbital Kondo effect can be realized in a quantum dot even in the 
absence of an a priori conserved quantum number like spin. The crossover 
between the interference- and Kondo-dominated regimes upon increase of 
the AB-phase <fi is most clearly seen in the single particle spectrum (Fig. ||) 
where the Kondo resonance gradually develops as <f> — > tt. Experimentally, 
the absorption power may be more easily accessible: With increasing </>, the 
absorption maximum is shifted from 5e to min(Tff , <5e), while at the same 
time the absorption intensity strongly increases (see Fig. |^). 




<o/r 



Fig. 8. Left: Effect of the AB-phase <j> on the single-particle spectrum. The total 
spectral density is shown for ei = — 1.6T, <S e = 0, U = 8. IT, and T = 0. Insets: 
Partial spectral densities for levels ci (left) and C2 (right). Same parameters as 
above, but with a finite level splitting Se = 0.08T. Right: Absorption power for 
U = 33T, ei = -3.3T, Se = 0.16T, T = and different values of <j>. 

So far, we have assumed that both levels are equally broadened by the 
reservoir, i.e. Ti = 2~2. In experiments, this need not be the case, although 
tuning within ST = 20% seems feasible. In order to clarify whether the orbital 
Kondo effect discussed previously is still visible under these conditions, we 
have therefore determined the influence of an asymmetric broadening on the 
single particle spectrum (Fig. ||). In the Kondo regime « n our model 
then corresponds to an effective Anderson Hamiltonian with spin-dependent 
hybridization, which is interesting in itself and has not been studied before. 
For the parameters used, we find a robust Kondo peak, which is split, but 
remains clearly visible even at an asymmetry of ST = 100%. We therefore 
conclude that the orbital Kondo effect discussed here should be accessible 
under realistic experimental conditions. 
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Fig. 9. Effect of an asymmetric broadening of the two levels for U = 16. 2A, ei = 
— 3.2-Ti, 4> = 7r, T = and different values of i~2. For comparison, the influence of 
a finite level splitting is also shown. 



Finally, we would like to extend our model to include the electron spin, 
which should give rise also to a magnetic Kondo effect. In this case, the dot 
Hamiltonian has to be generalized as 



-Hdot = 



+ Uriin 2 + U\m^n\i + U^n^n^ 



(14) 



where Uj = ^Za nia - F° r ^ ne s P ecia l case of U = U\ = U2 and Se = 0, 
this model has been studied before 27,28[], In realistic double-dot systems, 
however, we expect U < Ux,2- This is the parameter regime we will address 
here. 

NRG calculations for the model including spin are very expensive from 
a computational point of view. This is due to the larger size of the conduc- 
tion band Hilbert space, which now contains four instead of two additional 
fermionic degrees of freedom per iteration, as shown in Fig. h(J. In our calcu- 
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Fig. 10. Effective NRG Hamiltonian for the spin-degenerate dot (|l4|). 



lations, we kept about up to 1000 levels in each iteration step, which for 20 
iterations required about 10 hours of CPU time on an IBM Power 3. 
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Results are shown in Fig. Ill]: For vanishing AB-phase (f>, the two dot 
levels are coupled to the same reservoir. The resulting RKKY interaction 
[ p9| leads to an effective ferromagnetic exchange coupling between the levels. 
For zero interdot correlation (U = 0), the resulting side peaks dominate the 
single-particle spectrum (left plot). In addition, the screening of the total dot 
spin leads to the Kondo resonance at the Fermi level, which is suppressed for 
increasing broadening r. 
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Fig. 11. Partial spectral density for level ci of a spin-degenerate double quantum 
dot. Left: U = 0, Ui = U 2 = 2.0, ei = e 2 = -1.0 and <P = T = with varying T. 
Right: Same parameters, but now F — 0.61 is fixed and the interdot interaction U 
is tuned to different values. 



Switching on the interdot U as shown in the right plot of Fig. |ll| enhances 
the lower side peak. At U — 1, a discontinuous change in the ground state 
occupation number occurs, while at even larger values of U, one again obtains 
the interference "shoulder" already discussed in the spinless case. 

For a finite AB-phase (f> we expect interplay of the orbital and magnetic 
Kondo effects, which will be the subject of a future publication | p0[ . 

4 Conclusion 

In this paper we have discussed Wilson's Numerical Renormalization Group, 
a powerful nonperturbative method designed specifically to calculate low- 
temperature properties of quantum impurity systems. It is the only method 
that yields reliable results for systems with very different energy scales (small 
Kondo temperature Tk, large bandwidth). Recent experimental progress in 
the fabrication of ultrasmall quantum dots has made the preparation of ar- 
tificial Kondo "atoms" with well-controlled parameters possible. NRG is the 
method of choice for the theoretical interpretation of spectral and transport 
measurements in terms of single impurity models. 

We have presented an extended NRG algorithm (DM-NRG) suitable 
for calculating low-temperature dynamics in the full frequency range. This 
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method has been applied to calculate the spectrum of a quantum dot in a 
magnetic field. 

Furthermore, we have studied the interplay between interference and 
Kondo correlations in multi-level quantum dots. Orbital Kondo screening 
has been observed which can be tuned by an external Aharonov-Bohm phase 
and should be most easily visible in the optical absorption power. Additional 
spin degeneracy gives rise to the (Spin-) Kondo effect and causes an effec- 
tive RKKY-interaction between the two dot levels. We expect competition 
between magnetic and orbital screening in the presence of a finite AB-phase. 
These results should be useful for the interpretation of recent experiments on 
vertically coupled dots (see e.g. 31 ) or lateral multi-dot arrangements J32|. 

In future applications of the NRG method, more complex impurities will 
be considered - one may even try to model the lowest exitations of a com- 
plete molecule in order to describe recent transport experiments, see e.g. J33|. 
From the methodical point of view, the extension of NRG to non-equilibrium 
calculations remains a major challenge. 

The author would like to thank R. Bulla, H. Schoeller, and D. Vollhardt for 
valuable discussions. This work was supported by the Deutsche Forschungs- 
gemeinschaft through SFB 484. 
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